Accepted by ApJ on June 28, 2011 

Preprint typeset using WT^^ style emulatcapj v. 03/07/07 



O 

(N 

O 

u 

6 

> 

O 

rn 
O 



X 



EVOLUTIONARY TRACKS OF TIDALLY STIRRED DISKY DWARF GALAXIES 

EwA L. LoKAs\ Stelios Kazantzidis^ and Lucio Mayer^ 

Accepted by ApJ on June 28, 2011 

ABSTRACT 

Using collisionless iV-body simulations, we investigate the tidal evolution of late-type, rotationally 
supported dwarfs inside Milky Way-sized host galaxies. Our study focuses on a wide variety of dwarf 
orbital configurations and initial structures. During the evolution, the disky dwarfs undergo strong 
mass loss, the stellar disks are transformed into spheroids, and rotation is replaced by random motions 
of the stars. Thus, the late-type progenitors are transformed into early-type dwarfs as envisioned by 
the tidal stirring model for the formation of dwarf spheroidal (dSph) galaxies in the Local Group. We 
determine the photometric properties of the dwarfs, including the total visual magnitude, the half-light 
radius and the central surface brightness as they would be measured by an observer near the galactic 
center. Special emphasis is also placed on studying their kinematics and shapes. We demonstrate 
that the measured values are biased by a number of observational effects including the increasing 
angle of the observation cone near the orbital pericenter, the fact that away from the pericenter the 
tidal tails are typically oriented along the line of sight, and the fact that for most of the evolution 
the stellar components of the dwarfs are triaxial ellipsoids whose major axis tumbles with respect to 
the line of sight. Finally, we compare the measured properties of the simulated dwarfs to those of 
dwarf galaxies in the Local Group. The evolutionary tracks of the dwarfs in different parameter planes 
and the correlations between their different properties, especially the total magnitude and the surface 
brightness, strongly suggest that present-day dSph galaxies may have indeed formed from late-type 
progenitors as proposed by the tidal stirring scenario. 

Subject headings: galaxies: dwarf - galaxies: Local Group - galaxies: fundamental parameters - 
galaxies: kinematics and dynamics - galaxies: structure 



1. INTRODUCTION 

The population of dwarf galaxies of the Local Group 
(for a review see Mateo 1998; Tolstoy et al. 2009) can 
be divided into two rather distinct subgroups, the dwarf 
irregular (dirr) galaxies occupying larger distances from 
the big galaxies, the Milky Way and Andromeda, and 
dwarf spheroidal (dSph) objects clustering near their re- 
spective hosts. While the former are typically disky, 
bright, dominated by rotation, gas-rich and still forming 
stars, the latter are rounder, faint, supported mainly by 
random motions, gas-poor and dominated by old stellar 
populations. 

The dIrr galaxies are believed to be of primordial ori- 
gin, i.e. formed by gas accretion onto potential wells of 
small dark matter halos. The formation of dSph galax- 
ies, on the other hand, probably required some evolu- 
tion that may occur either in isolation by purely bary- 
onic processes such as cooling, star formation, feedback 
from supernovae and UV background radiation (Ricotti 
& Gnedin 2005; Tassis et al. 2008; Sawala et al. 2010) 
or be induced by gravitational interactions with the en- 
vironment in the form of galaxy harassment (e.g., Moore 
et al. 1996) resonant stripping (D'Onghia et al. 2009) 
or tidal stirring (Mayer et al. 2001a, b; Kazantzidis et al. 
2004). 
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The tidal stirring scenario proposed by Mayer et al. 
(2001a, b) is particularly promising as it provides a link 
between the two types of dwarf galaxies observed in the 
Local Group. It postulates that dirrs were actually pro- 
genitors of dSph galaxies and provides a mechanism for 
such transformation in the form of tidal stripping and 
tidal shocking of the dwarfs on orbits around the hosts. 
As the efficiency of these mechanisms increases with the 
decreasing distance of the dwarf galaxy from its host, it 
also naturally explains the morphology-distance relation 
observed in the Local Group. When combined with ram 
pressure stripping and the effect of the UV background 
(Bullock et al. 2000; Susa & Umemura 2004) the model 
can also account for the low gas content and the very 
high mass-to-light ratios of some dwarfs (Mayer et al. 
2006, 2007; see also Mayer 2010 for a review). 

The efficiency of the tidal transformation of disky 
dwarfs into dSph galaxies has been recently explored in 
great detail by Kazantzidis et al. (2011, hereafter Kll) 
for different orbital and structural configurations of the 
dwarfs interacting with the Milky Way. The dwarfs were 
initially composed of a stellar disk embedded in a more 
extended dark matter halo. They confirmed the general 
picture sketched earlier by Mayer et al. (2001a, b) and 
Klimentowski et al. (2009a): in most cases the tidal 
interaction of the dwarf with the Milky Way results in 
strong mass loss, the morphological transformation of the 
disk to a bar and then a spheroid and the transition from 
ordered (rotation) to random motion of the stars. Only 
in three out of 19 simulated cases (where the orbits were 
very extended or the dwarf's halo very concentrated) did 
the end product retain enough of its initial disky charac- 
teristics that it could not be classified as a dSph galaxy. 
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Similarly, in almost all cases, the morphological trans- 
formation involved the formation at the first pericenter 
passage of a tidally induced bar that was shortened in 
the subsequent evolution. 

The study of Kll focused on the intrinsic properties of 
the dwarfs which are best suited to describe their dynam- 
ics as precisely as possible, but can only be measured ac- 
curately in the idealized situation of A'^-body simulations 
where full 3D information is available. Here, we focus 
on strictly observational parameters that can be directly 
compared to real data. The purpose of the present work 
is therefore two-fold. First, we will measure the obser- 
vational parameters of the dwarfs in a realistic way and 
identify possible biases affecting the measurements in the 
environment characteristic for dwarf galaxies in the Local 
Group. Second, we will explore the correlations between 
the observed parameters of the dwarfs in order to see if 
these can provide some insight into their formation sce- 
narios. 

The paper is organized as follows. In section 2 we 
briefly describe the simulations used in this work. Sec- 
tion 3 discusses the way we measure the observational pa- 
rameters focusing on different observational biases that 
affect them. We also describe their variability on dif- 
ferent timescales and its origin. In section 4 we present 
long-term evolution of the measured parameters for dif- 
ferent orbits and different initial structure of the dwarfs. 
In section 5 we compare our results to observations. The 
discussion follows in the last section. 

2. THE SIMULATIONS 

In this study, we use a set of the collisionless A'^-body 
simulations presented in detail by Kll. The goal of the 
Kll study was to elucidate the formation of dSph galax- 
ies via tidal interactions between late-type, rotationally 
supported dwarfs and Milky Way-sized hosts. 

Kll employed the method of Widrow & Dubinski 
(2005) to construct numerical realizations of fully self- 
consistent dwarf galaxy models composed of exponen- 
tial stellar disks embedded in cuspy, cosmologically mo- 
tivated Navarro et al. (1996, hereafter NFW) dark mat- 
ter halos. The reference dwarf galaxy model had a virial 
mass of Mh = 10^ Mq and a concentration parameter 
c = 20. The disk mass fraction, nid, and the halo spin 
parameter. A, were equal to 0.02 and 0.04, respectively. 
The resulting disk radial scale length was i?d — 0.41 kpc 
(Mo et al. 1998) and the disk thickness was specified by 
the thickness parameter z^/Rd = 0.2, where denotes 
the (sech^) vertical scale height of the disk. We refer the 
reader to Kll for a detailed discussion regarding these 
choices. 

This default dwarf galaxy model was placed on seven 
different orbits around a single, massive host with the 
present-day structural properties of the Milky Way. In 
particular, the primary galaxy model was constructed as 
a live realization of the MWb model in Widrow & Dubin- 
ski (2005), which satisfies a broad range of observational 
constraints for the Galaxy and consists of an exponential 
stellar disk, a Hernquist (1990) bulge, and an NFW dark 
matter halo. The orbital apocenters, rapo, and pericen- 
ters, Tperi, of the dwarf around the host galaxy are listed 
in the third and fourth column of Table [l] and their 
choices are discussed thoroughly in Kll. The trajecto- 
ries of the reference dwarf galaxy model on the seven 
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Fig. 1. — Orbits of the dwarf galaxy in simulations 01-07. 
In all panels we repeat the orbit Ol for reference. In the upper 
panel all three orbits have the same eccentricity rapo/^peri = 5 
(typical for Milky Way subhalos) but different size. In the middle 
panel all three orbits have apocenter rapo = 125 kpc but different 
eccentricity. The lower panel shows the least (06) and the most 
(07) eccentric orbit in comparison to the typical one (Ol). Note 
that the scale of the middle panel is different than in of the other 
two. 
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TABLE 1 

Properties of the simulated dwarfs. 
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orbits are shown in Figure [T] projected onto the initial 
orbital plane. 

The simulations marked as 01-05 in Table [1] corre- 
spond to runs R1-R5 in Kll. Runs 06-07 are two addi- 
tional experiments mentioned in section 5.4 of Kll and 
discussed in more detail in Lokas et al. (2010c). Lastly, 
the alignments of the internal angular momentum of the 
dwarf, that of the primary disk and the orbital angular 
momentum were all mildly prograde and equal to z = 45° 
(see Kll for a discussion regarding the implications of 
this choice). Therefore, the results of the present study 
should not be affected by any strong coupling of angular 
momenta. Initial conditions for all numerical simulations 
were generated by building models of dwarf galaxies and 
placing them at the apocenters of their orbits. The tidal 
evolution of the dwarfs inside their host galaxy was fol- 
lowed for 10 Gyr using the multistepping, parallel, tree 
A^-body code PKDGRAV (Stadel 2001). 

In simulations S6-S17 (see Table [Ij, we varied the 
structural parameters of the dwarfs while keeping them 
on the same orbit (as in run 01 with rapo — 125 kpc and 
^'peri = 25 kpc). These experiments correspond to runs 
R6-R17 in Kll. The second column of Table [T] lists the 
parameters which were varied in each case; in parenthe- 
ses we give the value by which a given parameter was 
changed. Columns 5 and 6 of the Table list the orbital 
times and the time of the last apocenter for each simu- 
lation. 

Lastly, each dwarf galaxy model contained a total of 
2.2 million particles {N^ = 10^ dark matter particles and 
Nd = 1.2 X 10^ disk particles). The gravitational soften- 
ing was set to eh — 60 pc and ed — 15 pc for the particles 
in the two components, respectively. In addition, the 
simulations analyzed here use iVn = 10^ particles in the 
disk, iVe = 5 X 10^ in the bulge, and iVn 2 x 10^ in 
the dark matter halo of the host galaxy MWb, and em- 
ploy a gravitational softening of cd — 50 pc, eb — 50 pc, 
and ch = 2 kpc, respectively. The choice for the fairly 
large softening in the dark matter particles of the pri- 



mary galaxy was motivated by our desire to minimize 
discretness noise in the host potential. 

3. OBSERVING THE DWARFS 

3.1. The absolute magnitude 

In order to estimate the observational parameters of 
the dwarfs we place an imaginary observer at the center 
of the Milky Way and introduce a spherical coordinate 
system centered on this point, with the radial coordinate 
along the line connecting the center of the Milky Way 
with the dwarf and the angles corresponding to celes- 
tial coordinates usually applied in astronomy. Since the 
dwarfs typically lose a substantial fraction of their stars 
the measurement of the total magnitude of the dwarf 
galaxy as a function of time is non-trivial as it is not a 
priori obvious which stars should be included as mem- 
bers of the dwarf at a given instant of the evolution. In 
real observations this measurement is done by integrating 
light up to a given radius or up to infinity if an analytical 
approximation of the light distribution is used. Here we 
adopt a simple approach of counting all stars included 
within some fixed limiting radius rum. We adopt the 
value of this radius to be the distance from the center 
of the dwarf of the most distant star in the disk at the 
starting point of the evolution. The values of this radius 
are listed in the seventh column of Table [T] This choice 
is physically well justified: all the stars that are beyond 
this radius at any later stage were ejected due to tidal 
forces. 

Examples of the evolution of the total magnitude of the 
dwarf on orbits 01-07 measured within radii r < 
from the center of the dwarf as a function of time are 
illustrated by dashed lines in Figure [21 The values were 
obtained by counting the stars within rii„i, multiply- 
ing by the stellar mass and converting to the luminos- 
ity assuming a variable mass-to-light ratio for stars of 
M/Lv = (1 + O.15t)M0/L0 in the visual band, where 
t is the time from the start of the simulation, so that 
M/Lv — IM0/L0 at the start of the simulation (10 Gyr 
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Fig. 2. — The evolution of the total visual magnitude in time 
for simulations 01-07. The dashed line shows measurements in- 
side a sphere of radius centered on the dwarf where is the 
limiting radius containing all stars at the beginning of the simula- 
tion. The solid line shows measurements for an observer situated 
at the galactic center within the angular aperture corresponding 
to the projected radius R = r^^. Vertical dotted lines indicate 
pericenter passages. 



ago) and AI/Ly = 2.5M0/L0 at the end of the sim- 
ulation (which we assume to correspond to the present 
time). This choice approximates well the evolution of the 
simple low-metallicity stellar population in the standard 
model described by Bruzual & Chariot (2003, their Fig- 
ure 1). In this way we take into account the fading of 
the stellar population in time. 

The evolution of the total magnitude determined in 
this way is very regular. Until the first pericenter pas- 
sage in all cases the number of stars (and the stellar 
mass) remains constant and equal to the initial value 
(almost no stars are stripped beyond nim) so the de- 
creasing magnitudes are due entirely to the fading of the 
stellar population. In the runs on more eccentric orbits 
(01, 02, 04 and 07) soon after each pericenter passage 
(marked by vertical dotted lines in Figure [2]) the lumi- 
nosity decreases due to the ejection of stars beyond 
by tidal forces, then declines more slowly again until the 
next pericenter. For the other runs this dependence on 
the orbital position would only be seen if we adopted a 
constant AI /L for the stars while in Figure [5] it is domi- 
nated by the fading of the stellar population. Note that 
the vertical scale is different in each panel of the Figure, 
e.g. many more stars are lost on the tightest orbit (02) 
and the one with the smallest pericenter (04) compared 
to the reference orbit 01 while very few are ejected on 
the most extended orbit 03. 

An observer near the galactic center will obviously not 
be able to measure this intrinsic luminosity of the dwarfs 
but rather will count all the stars within the observation 
cone with the opening angle corresponding to at the 
distance of the dwarf. We assume that the distance of 
the dwarf is known perfectly (without any error) and 
the observer counts all stars within the projected radius 
R < ^"11111 corresponding to the angular distance from the 
center of the dwarf in the plane of the sky, which initially 
contained all the stars. The corresponding results for or- 
bits 01-07 are shown in Figure[2]as solid lines. For runs 
03, 05 and 06, where the fading dominates, the result 
is almost identical, but for the other cases this measure- 
ment is biased with respect to the intrinsic one described 
above. Both measurements coincide at the pericenters of 
the orbits, but in other parts of the orbit (except for the 
initial part of the evolution) the estimated magnitudes 
are brighter for the realistic observations. This is due to 
the fact that near pericenters the tidal tails are oriented 
along the orbit, i.e. perpendicular to the line of sight of 
the observer, while near apocenters their orientation is 
more along the observer's line of sight. This effect, de- 
scribed in more detail by Klimentowski et al. (2009b), is 
significantly stronger for orbits with higher eccentricity. 

The second observational effect is the presence of the 
spikes, sudden increases of the measured magnitude at 
the later pericenter passages, well visible especially for 
02 (second panel of Figure[2]). These are the result of the 
particularly strong stellar mass loss for this tight orbit 
combined with a large opening angle of the observation 
cone at pericenters (the distance of the dwarf from the 
galactic center is then about 17 kpc while its assumed 
diameter is 12.5 kpc). At later times there is a lot of 
debris around the Milky Way stripped from the dwarf at 
earlier passages (see the orbit plotted with the red line 
in the upper panel of Figure [ij which are included in 
the observer's field of view. Note that we include here 
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all the debris stars without discriminating by distance 
(since distances are typically not known for individual 
stars) and most often the older debris lies quite far from 
the dwarf itself. These stars, if far enough, could in prin- 
ciple be excluded from the sample by careful analysis of 
the color-magnitude diagrams of the real dwarfs as they 
would cause a spread in luminosity. This effect would 
however have to be disentangled from the effects of star 
formation history (see e.g. the discussion for Leo I in 
Gallart et al. 1999, Sohn et al. 2007 and Mateo et al. 
2008). 

Given that galaxies on eccentric orbits spend a small 
fraction of the orbital time at pericenters and such small 
pericenter distances are probably not common among 
Milky Way satellites (e.g. Lux et al. 2010) this effect 
should be rare. It could in principle be observed e.g. for 
the Sagittarius dwarf if it had completed many pericen- 
tric passages around the Milky Way, which is unlikely 
(see Lokas et al. 2010b). In the following we will use the 
measured values at apocenters, where dwarfs are most 
likely to be seen, as most representative for the observed 
magnitudes. At these parts of the orbits the differences 
between the values measured by our observer {R < rnm) 
and the intrinsic ones (r < rii,„) are never larger than 
0.4 mag (for run 04) so the tidal tails do not seem to 
bias the measurements very strongly even for the most 
heavily stripped dwarfs. 

3.2. The half-light radius 

The next structural parameter of special interest is the 
half-light radius ri/2. We first discuss the half-light ra- 
dius measured in 3D, i.e. the radius containing half the 
total luminosity calculated from the number of stars con- 
tained within the 3D fixed radius r < r^^. Half-light 
radii determined in this way are plotted as dashed lines in 
Figure |3] for orbits 01-07. The general trend, best seen 
for the tightest orbit 02, is that ri/2 decreases with time 
on a large timescale, which reflects the decreasing size of 
the dwarf due to tidal stripping. Significant temporary 
increase of rx/2 is however visible after each pericenter 
passage (indicated by vertical dotted lines in Figure [3]), 
especially for runs 04, 07 and 02 with the smallest peri- 
centers. This is due to the expansion of the dwarf's stel- 
lar component in response to the tidal forces which are 
strongest at the pericenter. 

A realistic measurement of ri/2, by an observer located 
at the center of the Milky Way, would be done in 2D, 
on the surface of the sky and using the total luminosity 
determined within R < rum. These measurements are 
shown in Figure |3] as solid lines. The half-light radii de- 
termined in this way are typically lower that their 3D 
counterparts which is due to the fact that the projected 
density profile of the stars is steeper, i.e. in the center 
of the dwarf all stars along the line of sight contribute 
and not only those that are close to the center in 3D. 
At later pericenters of simulation 02 the values of the 
projected rx/2 show strong peaks. This is a direct con- 
sequence of the peaks in the measured total magnitudes 
due to large opening angles of the observation cone at 
these instances, as discussed above. Another feature of 
the evolution of the projected ri/2 are wider bumps be- 
tween pericenters, especially well visible in run 04 (but 
also present in 02), where the measured values of pro- 
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Fig. 3. — The evolution of the half-Hght radius in time for simula- 
tions 01-07. The dashed line shows measurements inside a sphere 
of radius ri[^ centered on the dwarf where riim is the limiting radius 
containing all stars at the beginning of the simulation. The solid 
line shows measurements for an observer situated at the galactic 
center within the angular aperture corresponding to the projected 
radius R = rum. Vertical dotted lines indicate pericenter passages. 
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jected ri/2 are about twice larger at apocenter than at 
pericenter. These are due to increased number of stars 
from the tidal tails that for this highly eccentric orbit are 
oriented along the observer's line of sight for most of the 
time. 

Another interesting observational effect that manifests 
itself in the measured values of the 2D half-light radii 
is their strong variability or even oscillation at the early 
stages of the evolution, especially visible after the first 
and second pericenter. As discussed in detail by Kll, 
in most simulations studied here after the first pericen- 
ter the stellar disk transforms into a triaxial shape, usu- 
ally a prolate spheroid or a bar which, as the evolution 
proceeds, becomes more spherical. The tumbling of the 
ellipsoid causes the observed variation of the measured 
rxji: when the observer's line of sight is perpendicular 
to the major axis of the stellar component the dwarf ap- 
pears more extended and the measured half-light radius 
is larger; when the major axis is aligned with the ob- 
server's line of sight, the measured ri/2 is smaller. This 
effect is strongest for simulation 07, but present even in 
the case of run 03 where a bar does not form and the 
stellar component is oblate, but still triaxial. Note that 
in the case of simulations 02 and 04, where the dwarf 
is particularly strongly stirred and its stellar component 
becomes almost spherical after the third pericenter, this 
variability of is no longer seen at these later stages 
and the variation of ri/2 follows that of the total magni- 
tude. 

This interpretation of the measurements is confirmed 
by Figure 2] where we plot the angle between the major 
axis of the stellar component and the fixed x axis of the 
simulation box which lies in the orbital plane of all sim- 
ulations (at the initial configuration). The direction of 
the major axis was determined using all stars within a 
fixed radius of 0.6 kpc which corresponds to (1 — 2)ri/2 
in all cases. At all times we measured the angle between 
the X axis and the dwarf's major axis on the side of the 
dwarf which is closer to the x axis so all angle values are 
between 0° and 90°. The values of the angle were plotted 
starting from the first pericenter when the shape of the 
dwarf's stellar component changes from the disk to the 
triaxial ellipsoid so that the major axis is well defined. 
For orbits 02 and 04 we do not show the measurements 
until the end because on these tight orbits the dwarf be- 
comes almost spherical early on and the major axis is 
then no longer well defined again. The amplitude of the 
angle changes in some cases, which is due to the preces- 
sion of the rotation axis. In addition, the rate at which 
the angle changes increases after the second pericenter 
for 01 which is due to speeding up of the figure rotation 
of the ellipsoid in this case, caused by the tidal force, as 
discussed by Kll (their section 5.2). 

3.3. The central surface brightness 

Another characteristic photometric property of dwarf 
galaxies is the central surface brightness fj,v- This pa- 
rameter does not have any 3D analog so we measure it 
directly in 2D by counting stars within O.lri/2 (where 
ri/2 is the half-light radius determined in 2D) and divid- 
ing by the surface of the circle of this radius. This scale 
is about the closest to the center of the dwarf possible, 
i.e. the position of the innermost data point in measure- 
ments of surface brightness in real data. The results of 
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Fig. 4. — The evolution of the angle between the major axis 
of the stellar component of the dwarf and the fixed x axis of the 
simulation box in time for simulations 01-07. Measurements are 
performed for stars within a fixed radius r < 0.6 kpc. Vertical 
dotted lines indicate pericenter passages. 

this procedure are shown as solid lines in Figure [5] for 
runs 01-07. 

We notice a similar degree of variability in the mea- 
sured surface brightness values as was observed for the 
2D half-light radii. Since we measure the surface bright- 
ness within a fixed fraction of ri /2 the two are obviously 
tightly related. Note that in the case of simulations 02 
and 04 the strong variability at the early stages is re- 
placed by an almost constant value between pericenters 
later on when the dwarf becomes spherical. The general 
trend is however to decrease the central surface bright- 
ness which means that the stars are stripped from all 
radii, not only from the outer parts i.e. the density pro- 
file of the stars is shifted down at all scales due to tidal 
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Fig. 5. — The evolution of the central surface brightness in 
time for simulations 01-07. Measurements are performed within 
O.lrj/2 where determined by an observer situated at the 

galactic center within the angular aperture corresponding to the 
projected radius R = r^^. Vertical dotted lines indicate pericenter 
passages. 



stripping. 



3.4. Kinematics 



One of the key parameters that distinguish dirrs from 
dSphs is their V/a value, where V is the rotation veloc- 
ity and cr is the (central) velocity dispersion of the stars. 
While dIrrs are believed to be rotationally supported, 
with V/a of the order of a few, dSphs usually exhibit 
much lower rotation levels with V/a below unity. How- 
ever, contrary to common belief, dSphs are not in general 
completely devoid of rotation and it has been detected for 
many dSph galaxies, including Ursa Minor (Hargreaves 
et al. 1994; Armandroff et al. 1995), Carina (Mufioz et 




Fig. 6. — The distribution of stars in simulation Ol within 2r-i/2 
as seen from the galactic center at the last apocenter of the orbit 
(green dots). The position of the stars were rotated so that the 
major axis of the galaxy image is along the x axis. Only 1 percent 
of the stars is shown for clarity. The inner black circle of radius 

0. 1rj/2 encloses stars that are used in the determination of the 
central velocity dispersion. Black rectangles indicate eight bins 
along the major axis of the dwarf in which the rotation velocity 
and the velocity dispersion profiles are measured. 

al. 2006), Sculptor (Battaglia et al. 2008), Leo I (Sohn 
et al. 2007; Mateo et al. 2008), Cetus (Lewis et al. 2007) 
and Tucana (Fraternali et al. 2009). 

We measured V and a using the following procedure, 
as close as possible to real observations. For each sim- 
ulation output we selected stars, as seen by an observer 
at the galactic center, within projected radii R < 2ri/2 
from the center of the dwarf, where /2 are the 2D half- 
light radii determined above. This is approximately the 
region where kinematic measurements for dwarf galaxies 
are typically performed. We then determined the prin- 
cipal axes of the 2D distribution of the stars and ro- 
tated their positions in order to align the major axis of 
the dwarf image with the horizontal (x) axis. Figure [6] 
shows the distribution of (1 percent) of the stars within 
R < 2ri/2 obtained in this way for the configuration cor- 
responding to the last apocenter for simulation 01. 

In order to estimate the central velocity dispersion of 
the dwarfs we selected for each output the stars within 
R < O.lri/2 (small black circle in Figure 16]) and cal- 
culated the dispersion of radial velocity using standard 
estimators (see e.g. Lokas et al. 2005). By the choice of 
this innermost region we make sure that the sample is not 
contaminated by tidally stripped stars in the tails. Next, 
we select a narrow strip of stars along the major axis of 
the dwarf of size O.lri/2 ^ 2ri/2 and sample the velocities 
of the stars in eight bins of size O.lri/2 x 0.5ri/2 each. 
Since the outer bins are likely to be contaminated by 
tidally stripped stars we introduce a cut-off in velocities, 

1. e. we exclude the stars with velocities that differ from 
the dwarf's systemic velocity by more than 3cr where a 
is the central velocity dispersion measured in the inner 
R < O.lri/2- Such a conservative cut-off is usually not 
sufficient for the reliable selection of stars for dynamical 
modelling (see Klimentowski et al. 2007, 2009b; Lokas 
et al. 2008; Lokas 2009) but seems stringent enough to 
measure basic kinematic properties of the dwarfs, as we 
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Fig. 7. — The velocity dispersion profiles (solid lines), the rota- 
tion velocity profiles (dashed lines) and the ratios of the two (dotted 
lines) measured along the major axis of the observed dwrarf for runs 
01-07 at their 2nd (left column) and last (right column) apocen- 
ters (except for the case 03 where the 2nd apocenter is the last so 
we plot the result for the final simulation output). 



aim to do here. 

The results of these measurements for simulations 01- 
07 and S6-S17 are shown in Figures [7] and |S] respectively 
for the outputs corresponding to the 2nd and the last 
apocenter (the left and right panel for each run). With 
dashed lines we plot the rotation velocity profiles and 
with solid lines the velocity dispersion profiles to illus- 
trate their dependence on the distance from the center 
of the dwarf. Dotted lines give the ratio of the two val- 
ues as a function of radius. All measurements were sym- 
metrized, i.e. we show the average of two measurements 
on both sides of the dwarf. In most cases, the rotation 
profile decreases significantly between the 2nd and the 
last apocenter, while the dispersion remains roughly on 
the same level. This is due to the fact that tidal stirring 
causes the rotation to be replaced by random motions 
increasing the velocity dispersion, but at the same time 
dwarfs lose mass and the dispersion is decreased. The 
net outcome of these two effects is that the dispersion 
remains roughly constant in time. 

All dispersion profiles look well-behaved, declining 
with radius, except for the case 04 (the right panel of 
the fourth row in Figure [71 orange line) where the disper- 
sion profile shows a secondary increase at larger distances 
from the center. This is due to the contamination from 
the tidal tails that affects our measurements because the 
nominal half-light radius determined for this case is very 
large. As already discussed above, and seen in the fourth 
panel of Figure O in this case the ri/2 values measured 
near apocenters at the later stages of evolution are about 
a factor of two larger than those at pericenters. Again, 
this is caused by the presence of tidal tails that for this 
very radial orbit (r'apo/^'pori = 10, see Table [1]) for most 
of the time (except at the very pericenter) are oriented 
along the observer's line of sight. The dwarf thus appears 
brighter and more extended (see also Figure 11 in Kli- 
mentowski et al. 2007). The contamination is also seen 
in the rotation curve (dashed orange line in the same 
panel of Figure [7]) which shows tidally induced velocity 
gradient in the outer bins while no rotation is detected 
in the inner part for this very evolved dwarf. 

Figure [9] shows the evolution of the kinematic prop- 
erties in time for the seven simulations 01-07 with dif- 
ferent orbits. In each panel the measurements of the 
central velocity dispersion are plotted with solid lines 
while the values of the rotation velocity are shown with 
dashed lines. The values of rotation are taken as a max- 
imum rotation velocity found along the rotation velocity 
profile (see Figures [7] and [8]) measured along the major 
axis within R < 2ri/2- In simulations 02 and 04 where 
the dwarfs are most strongly affected by tidal forces, the 
contamination by tidal tails is well visible, although in a 
slightly different way in each run. For orbit 02 we see 
sharp increases of both velocity dispersion and rotation 
velocity at later pericenters caused by spikes in the mea- 
sured values of ri/2. In simulation 04, on the other hand, 
while the measured values of the central velocity disper- 
sion remain unaffected (until the very last pericenter), 
the rotation velocity is overestimated when the dwarf is 
on its way from the pericenter to the apocenter. As dis- 
cussed above, this is due to the overestimated value of 
ri/2 which makes us probe the region well beyond the 
dwarf's main stellar body. 

Finally, it is worth noting that the tidal effects on the 
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Fig. 8. — The velocity dispersion profiles (solid lines), the rotation velocity profiles (dashed lines) and the ratios of the tv^o (dotted lines) 
measured along the major axis of the observed dwarf for runs S6-S17 at their 2nd and last apocenters. 



observables in runs 02 and 04, although in both cases 
due to tidally stripped stars, are of slightly different ori- 
gin; while in 02 they are due to earlier wraps of tidally 
lost material (this run has more pericentric passages than 
any other orbit we considered) , in 04 they are due to re- 
cently lost material residing in the tails oriented radially 
because of the high eccentricity of the orbit. In other runs 
the most characteristic feature of the kinematic measure- 
ments is the oscillation of the values (especially those 
of the central velocity dispersion) in simulations due, as 



discussed above, to the non-sphericity of the stellar dis- 
tribution and the tumbling of its triaxial shape. 

3.5. Shapes 

We conclude the measurements of the observational 
parameters of our simulated dwarfs by determining their 
shapes. For this purpose we select in each output the 
stars within projected radii R < 2ri/2 and rotate the 
2D distribution to align the major axis of the galaxy 
image with x and the minor axis with y, exactly as for 
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Fig. 9. — The evolution of the central velocity dispersion (solid 
line) and the rotation velocity (dashed line) in time for simulations 
01-07. The central velocity dispersion was measured for stars with 
projected radii R < ^-^id the rotation velocity was taken as 

a maximum value found along the major axis of the galaxy image 
within 2ri/2i shown in the upper panel of Figure |6] Vertical 
dotted lines indicate pericenter passages. 
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Fig. 10. — The evolution of the ellipticity e = 1 — b/a in time 
for simulations 01-07. Measurements of the axis ratios b/a are 
performed within 2ri/2 where r^/g is determined by an observer 
situated at the galactic center within the angular aperture corre- 
sponding to the projected radius R = rum- Vertical dotted lines 
indicate pericenter passages. 

the measurements of the kinematic properties. We then 
calculate the ellipticity parameter e — 1 ~ b/a where b/a 
is the axis ratio of the 2D distribution determined from 
the eigenvalues of the 2D inertia tensor. 

The evolution of the ellipticity parameter e in time is 
shown in Figure [10] for simulations 01-07. Note that 
all runs start from the same value of e = 0.33 which is 
due to the specific orientation of the dwarf galaxy disk 
at the beginning of the simulation. Namely, the disks 
are inclined by 45° to the initial orbital plane and the 
initial velocity vector so the galaxy images appear el- 
liptical to the observer at the center of the Milky Way. 
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Fig. 11. — The evolution of the total visual magnitude in time 
from the first to the last apocentcr for different orbits (simulations 
01-07, upper panel) and different dwarf models (simulations S6- 
S17, lower panel). 
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Fig. 12. — The evolution of the half light radius in time from the 
first to the last apocenter for different orbits (simulations 01-07, 
upper panel) and different dwarf models (simulations S6-S17, lower 
panel). 

At later stages the ellipticity parameter displays strong 
variability due to the formation of the triaxial shape and 
tumbling of this shape, as well as the variability of the 
2ri/2 scale inside which we measure the shape. A clear 
transition to sphericity is however seen for runs 02 and 
04 where the ellipticity is very close to zero at the later 
stages of the evolution with momentary exceptions at 
pericenters where the scale 2ri^2 picks up some of the 
tidal tails oriented perpendicular to the line of sight at 
these instants. 

4. LONG-TERM EVOLUTION OF THE OBSERVED 
PROPERTIES 

In order to study the general trends in the evolution of 
the observational parameters in time in Figures [TTII15I we 
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Fig. 13. — The evolution of the central surface brightness in time 
from the first to the last apocenter for different orbits (simulations 
01-07, upper panel) and different dwarf models (simulations S6- 
817, lower panel). 
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Fig. 14. — The evolution of the ratio V/cr in time from the 
first to the last apocenter for different orbits (simulations 01-07, 
upper panel) and different dwarf models (simulations S6-S17, lower 
panel). 

plot the values measured previously only at subsequent 
apocenters, but for all simulations. In the upper panel of 
each figure we show with different colors the results for 
the simulations for different orbits and the same initial 
dwarf structure 01-07, while in the lower panel we plot 
results for different initial structures of the dwarf and for 
the same default orbit S6-S17. A particular simulation 
can be identified by referring to the colors listed in the 
last column of Table [TJ The time of the last apocen- 
ter for each simulation can be found in the sixth column 
of the Table and the corresponding 2D values of the to- 
tal magnitude My, half-light radius ri/2, central surface 
brightness /iy, the ratio of the rotation velocity to the 
central velocity dispersion V/a and ellipticity e = 1 — b/a 
measured at this last apocenter are listed in columns 8-12 
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Fig. 15. — The evolution of the ellipticity e = l — b/a'm time from 
the first to the last apocenter for different orbits (simulations Ol- 
07, upper panel) and different dwarf models (simulations S6-S17, 
lower panel). 

of the Table. 

The general trends in the evolution of the parameters 
are now clearly visible. As expected, the tidal stripping 
decreases the total magnitude of the dwarf and most so 
for the tightest orbit 02 and the one with the smallest 
pericenter and reasonably small orbital time 04. The 
trend of decreasing with time is less clear and more 
visible for dwarfs of different structure rather than on 
different orbits. Overall, the half-light radius remains 
remarkably constant over time, in contrast to the other 
characteristic scale- length of the dwarf, rmax, where the 
maximum of the circular velocity occurs which was found 
by Kll to decrease very strongly during the evolution. 
This is not surprising, however, given that Tmax is a char- 
acteristic scale- length of the total gravitational potential, 
which is dominated by dark matter, while ri/2 charac- 
terizes the distribution of stars which are not so strongly 
affected by tidal forces. 

The central surface brightness behaves differently from 
the other two properties discussed above. Typically, its 
value increases between the first and the second apocen- 
ter and then starts to decrease. This is due to the for- 
mation of the bar which in most cases occurs at the first 
pericenter and increases the density of the stars in the 
center of the dwarf. Later on, when the stellar distribu- 
tion becomes more spherical and more stars are stripped, 
fiv decreases as a function of time. 

In Figure [T3] we plot the ratio of the rotation velocity 
and the central velocity dispersion, V/a, which measures 
the amount of ordered versus random motion of the stars. 
Both for different orbits and for different dwarf models 
the trend of V/ a decreasing with time is clearly visible. 
This is one of the key signatures of tidal stirring and con- 
firms that the process may indeed transform rotationally 
supported systems to those supported by random mo- 
tions. 

Finally, Figure [15] shows the trends in the evolution 
of the observed ellipticity. Here the trends are diffi- 
cult to detect since a similar small ellipticity is mea- 



sured for a disk close to face-on and a genuine almost 
spherical galaxy. In almost all our simulations the disk 
was initially inclined so that the measured ellipticity was 
e = 0.33. The only exceptions are runs S6 where the 
disk is seen edge-on resulting in a large initial elliptic- 
ity e = 0.81 (green line in the lower panel of Figure ITSj) 
and S7 where the disk is more face-on than the default 
orientation resulting in a small initial e = 0.18 (red line 
in the lower panel of Figure [T5|) . For different orbits the 
transition to the almost perfectly spherical shape is how- 
ever reflected in very low e values for runs 02 and 04, 
while for different dwarf models the ellipticity averaged 
over all runs is smaller than initially. 

5. COMPARISON WITH OBSERVATIONS 

The ultimate purpose of our study is to compare 
the measured observational parameters of the simulated 
dwarfs to real data. In Figures [TBllTQl we plot the central 
surface brightness ^y, the 2D half-light radius ri/2, the 
kinematic parameter V/a and the ellipticity e = 1 — b/a 
respectively as a function of the total magnitude of the 
simulated dwarfs at subsequent apocenters. The color 
coding of the simulations is the same as in previous Fig- 
ures (see the last column of Tabled]) and again the upper 
panel is for dwarfs on different orbits while the lower one 
for dwarfs of different structure on the same default orbit. 
The direction of the lines from the left to the right, i.e. 
from the brighter to the fainter magnitudes, corresponds 
to the time flow because the total magnitude becomes 
fainter with time in all cases. 

In Table [2] we list the properties of dwarf galax- 
ies of the Local Group with magnitudes in the range 
— 16 < My < —8 currently available in the literature. 
The columns give the galaxy name, the morphological 
type, the total visual magnitude, the half-light radius, 
the central surface brightness, the kinematic parameter 
V/a and the ellipticity parameter e = 1 — b/a. The 
references are provided in the last column of the Table. 
In cases where the exponential scale-lengths were only 
available, they were multiplied by a factor of 1.7 to ap- 
proximate the half-light radius of the exponential profile. 
For the V/a, the 'raw' values are used, i.e. those mea- 
sured directly from the data, as we do for our simulated 
dwarfs, without correcting for inclination. 

The circles in Figures [TBIIT^ show the real data for 
dwarf galaxies in the Local Group from Table [2] The 
open circles mark the data for dirr galaxies, filled black 
circles for the dSph and dSph/dE galaxies and gray cir- 
cles for transitory dIrr/dSph dwarfs. The data in the /xy- 
Mv and ri/2-Mv plots show strong correlation between 
the parameters (although in the latter case the scatter is 
much larger) demonstrating the trends of fainter galax- 
ies to possess lower surface brightness and smaller size. 
There is also an obvious trend in V/a- My- rotationally 
supported systems are typically brighter. The correla- 
tion is weakest in the e-My plane due to the difficulty 
in relating the observed ellipticity to the real shape of a 
galaxy, as discussed above. 

The evolutionary tracks of the simulated dwarfs in the 
Hv-Mv plane (Figure [TBI) have a characteristic shape, 
similar for all simulations. While the magnitudes evolve 
towards the fainter end monotonically, the values of the 
surface brightness between the first and second apocen- 
ter increase and then mildly decrease at all subsequent 
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TABLE 2 

Properties of the dwarf galaxies of the Local Group. 



Dwarf 


Type 


My 


1/2 




V/<j 


e = 1 — b/a 


References 


galaxy 




[mag] 


[kpc] 


[mag arcsec"'^] 








WLM 


"dlrr 


—14.5 


1.509 


23.6 


2.63 


0.59 


1 3 


IC 10 


dirr 


—15.7 


0.775 


22.1 


3.75 


0.3 


I 


NGC 147 


dSph/dE 


—15.5 


0.753 


21.6 


1.0 


0.46 


1 2 


NGC 185 


dSnh /dR 


—15.5 


0.521 


20.1 


0.63 


0.26 


1 2 


LGS 3 


dIrr/dSph 


-10.5 


0.312 


25.6 


0.22 


0.26 


1,3 


IC 1613 




— 14.7 


1.870 


23.5 


2.47 


0.24 


1,3 


Phoenix 


dlrr /dSnh 


— 10.1 






0.22 


0.3 


1 


UGCA 92 


dlrr 


-12.6 


0.441 


24.2 


4.13 


0.55 


1,4 


Leo A 


dlrr 


— 11.4 


0.341 




0.32 


0.40 


1,19 


Ssxtans B 


dlrr 


— 14.2 


0.417 


22.9 


1.22 


0.23 


13 4 


NGC 31 OQ 




— 15.7 


1.920 


23.6 


6.7 


0.80 


1 


Antlia 


dlrr/dSnh 


—10.8 


0.672 


24.3 




0.35 


1 






—14.6 


1.068 


23.7 


2.38 


0.21 


1 3 


GR 8 


dlrr 


—11.6 


0.189 


22.3 


0.64 


0.31 


1 
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Fig. 16. — The central surface brightness as a function of the 
total magnitude. Colored lines show the results for the simulated 
dwarfs at subsequent apocenters for different orbits (simulations 
01-07, upper panel) and different dwarf models (simulations SC- 
SI?, lower panel). The evolution proceeds from the left to the 
right, from the brighter to the fainter magnitudes in all cases. In 
both panels the open circles show the data for dlrr galaxies, filled 
black circles the data for the dSph and dSph/ dE galaxies and gray 
circles for transitory dlrr/dSph dwarfs. All data for the real dwarfs 
are listed in Table 1^ 

apocenters. This later evolution traces very well the cor- 
relation between the observed values of the parameters 
for real Local Group dwarfs and is especially well visible 
in the case of runs 02 and 04 (the red and orange line 
in Figure I16p where the dwarfs are most affected by tides 
and therefore most evolved. 

In the ri /2-Mv plane (Figure [T7| the evolution of the 
simulated dwarfs shows more variety. Typically, the 
characteristic radii decrease strongly between the first 
and second pericenter, but in the later evolution they 
can also increase. Interestingly, because of the contam- 
ination by tidal tails, it can also happen for the most 
affected dwarfs, like the one of run 04 (the orange line 
in the upper panel of Figure [T71) . where ri/2 increases 
during the later evolution, in contrast to the case of 02 
(the red line). Overall, different orbits and different ini- 



FlG. 17. — Same as Figure [161 but for the half-light radius as a 
function of the total magnitude. 

tial structures of the dwarf are able to reproduce the 
trends seen in the real data. 

The direction of the evolution of the simulated dwarfs 
in the ri/2-Mv and especially in the /iy-My plane 
strongly suggests that the observed correlations may be 
easily explained if dSph galaxies of the Local Group 
evolved from late-type dwarfs as predicted by the tidal 
stirring scenario. The trend of surface brightness de- 
creasing with decreasing luminosity as observed in real 
dSph galaxies, in contrast to ellipticals, has been inter- 
preted as pointing towards different formation scenar- 
ios of spheroidal and elliptical galaxies by Kormendy 
(1985) and Kormendy et al. (2009). They suggested 
that elliptical galaxies may form mostly via mergers while 
spheroidals are rather late-type systems that underwent 
transformation due to star-formation processes or envi- 
ronmental effects such as tidal stirring. This idea is sup- 
ported by the analysis of the simulated Local Group by 
Klimentowski et al. (2010) who found that mergers of 
subhalos are rather rare in such systems and occur early 
on so they cannot significantly contribute to the forma- 
tion of a large fraction of dSph galaxies. The tidal stir- 
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Fig. 18. — Same as Figure [161 but for the kinematic parameter 
y/iT as a function of the total magnitude. 

ring scenario thus seems to be the most effective gravi- 
tational mechanism by which such objects could form. 

6. DISCUSSION 
6.1. Observed versus intrinsic properties 

Our results demonstrate the difficulty in determining 
the intrinsic properties of tidally stirred dwarf galax- 
ies from the photometric and kinematic measurements 
alone. The interpretation of the observational biases 
in the photometric properties such as total magnitudes, 
half-light radii and central surface brightness is straight- 
forward. All of them are affected by three issues: 
the presence of contamination by tidal tails, the non- 
sphericity of the stellar distribution and the 'proximity 
effect' that for strongly evolved dwarfs on orbits with 
small pericenters makes the observer seeing the dwarf at 
pericenter include more stars than are actually associated 
with the dwarf. 

The measurements of kinematics and shapes are more 
complicated to interpret. Although in both these quanti- 
ties the trend towards smaller V/a (signifying less rota- 
tion and more random motion in the stars) and towards 
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Fig. 19. — Same as Figure [161 but for the ellipticity parameter 
e = 1 — 6/a as a function of the total magnitude. 

more spherical shapes is clearly seen especially for more 
evolved dwarfs, these measurements depend critically on 
the initial orientation of the dwarf disks. This orienta- 
tion was the same for all our simulations except S6 and 
S7 and may be considered 'typical' in the sense that our 
observer at the galactic center is able to detect both ro- 
tation and non-spherical shape of the stellar component 
at the initial state and later. 

In particular, the cases 03, 05 and S15 classified as 
'non-dSph' by Kll (see cases R3, R5 and R15 in their 
Table 2) based on their high intrinsic V/a > 1 (mea- 
sured as the mean rotation around the shortest axis and 
the mean velocity dispersion for stars inside rmax) also 
show high values of the order of unity in our measure- 
ments so they would be classified as transitory objects 
also by our observer. However, if the disk was seen ex- 
actly edge-on (as in run S6) the observer could identify it 
as rotating disk based on shape and rotation alone, while 
if it was face-on our observer would not be able to dis- 
tinguish the disk from a spheroid supported by random 
motions. In real observations such cases can be morpho- 
logically classified by looking for additional signatures 
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Fig. 20. — The central surface brightness (upper left panel), the half-light radius (lower left panel), the V/cr parameter (upper right 
panel) and the cUipticity (lower right panel) as a function of the total magnitude. Green squares show the measured parameters of the final 
products of tidal stirring in our simulations (at the last apocenter of their orbit). The open black circles show the data for dirr galaxies, 
filled black circles the data for the dSph and dSph/dE galaxies and gray circles for transitory dIrr/dSph dwarfs from Table [2] 



like the presence of gas or star formation regions. 

Despite these difficulties, we emphasize that in this 
study we measured the properties of the simulated dwarfs 
in a realistic way so that they can be directly compared to 
real data. Figure [501 summarizes our results by showing 
as green squares the measured properties of the simulated 
dwarfs at the last apocenter in the fj^y-My, f 1/2- My, 
Vja-My and e-My planes. The circles display the data 
for dwarf galaxies of the Local Group from Table [5] as in 
Figures 11611191 The comparison proves that our dwarfs 
are indeed very similar to the real ones in those crucial 
parameters and thus shows that tidal stirring indeed can 
produce dSph galaxies from dirr-like progenitors. Note 
that all our dwarfs (except two) initially had rather high 
masses of 10^ M0 and therefore ended up on the more 
luminous side of the observed distribution. Starting from 
lower masses and luminosities we could easily reproduce 
also the less luminous dwarfs. 

6.2. Comparison with earlier work 

Observational parameters of the dSph galaxies pro- 
duced by tidal stirring have already been discussed in 
the papers proposing the tidal stirring as a possible sce- 



nario for the formation of dSph galaxies by Mayer et al. 
(2001a, b). They have shown that the properties of the 
final products indeed follow the ^y — My relation found 
for real dwarfs and that the V/ a and ellipticity both de- 
crease with time. Lokas et al. (2010a) measured the ob- 
servational properties of dSph galaxies formed in three 
simulations discussed in Klimentowski et al. (2009a), 
similar to those studied here, and found them to depend 
strongly on the line of sight but still in the range char- 
acteristic of real dwarfs. In this work we improved on 
these earlier results by showing the actual evolutionary 
tracks of the dwarfs in different parameter planes and 
considering different observational biases caused by the 
observer's position close to the galactic center. 

Peharrubia et al. (2008) studied the effect of tides 
on objects that were spherical from the start and con- 
structed as stellar King profiles embedded in more ex- 
tended dark halos. Understandably, their dwarfs do not 
undergo any morphological transformation and remain 
spherical and supported by random motions with no net 
rotation. Thus, they were unable to address such issues 
as the evolution of rotation or ellipticity. Their analysis 
was therefore restricted to such observational parame- 
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ters as the total luminosity, the half-light radius, sur- 
face brightness and velocity dispersion. The evolutionary 
tracks of their spherical dwarfs in the /iy-My and rxji- 
My planes (their Figure 10) are similar as in our simu- 
lations (Figures ITBlfTT)) . i.e. in both cases there is a clear 
trend of surface brightness and radii decreasing with de- 
creasing luminosity of the dwarfs in time. Interestingly, 
the characteristic scales of the stellar distribution, quan- 
tified by the core radius of the King profile, were found 
there to evolve rather mildly and decrease only by about 
a factor of two, as in the case of our half-light radius. 

Despite these similarities in the photometric proper- 
ties, we find very different results concerning the evo- 
lution of the velocity dispersion and the mass-to-light 
[M / L) ratio. Since we were more interested in the evo- 
lution of the kinematic parameter V/ a we do not show 
analogous plots for the dispersion alone, but the corre- 
sponding trends can be read from Figure ^ where both 
quantities, V and cr, are plotted separately. Clearly, the 
central velocity dispersion decreases systematically only 
for the most strongly evolved dwarfs (02 and 04), while 
it remains roughly constant for the other cases, in con- 
trast to the findings of Peharrubia et al. The evolution of 
the M / L ratio as an intrinsic, rather than observational, 
parameter of our dwarfs was discussed already by Kll: 
they found that tidal stripping typically decreases M/L 
except for strongly evolved dwarfs where the dark halo 
has been truncated down to the stellar scales. Again, 
our results on this point are discrepant with those of 
Peharrubia et al. who found M/L increasing, making 
the dwarfs darker with time. The differences in the evo- 
lution are due to different initial conditions: in our case 
the stars are more tightly bound (especially after the 
formation of the bar) and are therefore more difficult to 
strip, keeping the velocity dispersion roughly constant 
and M/L rather low. In the King models of Peharrubia 
et al. the stars are weakly bound and easily stripped 
which decreases their dispersion and increases M/L. In 
summary, while the evolution of the photometric prop- 
erties seems to be a general feature of tidal stripping, 
largely independent of the initial structure of the dwarfs, 
the kinematic and dynamical properties depend on the 
assumed dwarf model. 

6.3. Conclusions 

We have studied the observational parameters of sim- 
ulated tidally stirred dwarf galaxies initially composed 
of disks embedded in extended dark matter halos. Our 
results build on and extend the previous study by Kll, 
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where intrinsic properties of the simulated dwarfs were 
discussed. Our main conclusions may be summarized as 
follows: 

1. The photometric parameters, such as the absolute 
magnitude, the half-light radius and the central 
surface brightness as well as the kinematic proper- 
ties can be reliably measured but are subject to a 
number of observational biases caused by the pres- 
ence of tidally stripped stars and the non-sphericity 
of the stellar component. 

2. The amount of bias in the measured quantities de- 
pends on the orbit and structural parameters of the 
dwarf but typically does not exceed 0.4 mag in total 
visual magnitude and a factor of two in half-light 
radius. Rotation can be overestimated for strongly 
stripped dwarfs on eccentric orbits even if a 3ct cut- 
off is applied to radial velocities. 

3. The effects of tidal stirring manifest themselves by 
decreasing the observed total magnitude, the half- 
light radius, the central surface brightness, rota- 
tion and ellipticity in time, while the central ve- 
locity dispersion remains roughly constant except 
for the most heavily stripped dwarfs where it also 
decreases due to strong mass loss. 

4. The correlations between the observational param- 
eters of the simulated dwarfs reveal strong simi- 
larity to those in the real data. The evolutionary 
tracks of the tidally stirred dwarfs move them from 
regions characteristic for dirr galaxies to those typ- 
ical for dSph galaxies of the Local Group. This be- 
havior corroborates the existence of a connection 
between these two types of objects and therefore 
strongly supports the tidal stirring model as a pos- 
sible scenario of their evolution. 
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